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ABSTRACT 

We study the physical properties derived from interstellar cloud complexes 
having a fractal structure. We first generate fractal clouds with a given fractal 
dimension and associate each clump with a maximum in the resulting density 
held. Then, we discuss the effect that different criteria for clump selection has on 
the derived global properties. We calculate the masses, sizes and average densities 
of the clumps as a function of the fractal dimension (Df) and the fraction of the 
total mass in the form of clumps (e). In general, clump mass does not fulfill a 
simple power law with size of the type M c i a W cl , instead the power changes, 
from 7 ~ 3 at small sizes to 7 < 3 at larger sizes. The number of clumps 
per logarithmic mass interval can be fitted to a power law N c i oc MJ“ M in 
the range of relatively large masses, and the corresponding size distribution is 
N c i oc at large sizes. When all the mass is forming clumps (e = 1) we 

obtain that as Df increases from 2 to 3 increases from ~ 0.3 to ~ 0.6 and 
an increases from ~ 1.0 to ~ 2.1. Comparison with observations suggests that 
Df ~ 2.6 is roughly consistent with the average properties of the ISM. On the 
other hand, as the fraction of mass in clumps decreases (e < 1) % increases and 
an decreases. When only ~ 10% of the complex mass is in the form of dense 
clumps we obtain ~ 1.2 for Df = 2.6 (not very different from the Salpeter 
value 1.35), suggesting this a likely link between the stellar initial mass function 
and the internal structure of molecular cloud complexes. 

Subject headings: ISM: structure — ISM: clouds — ISM: general 


1. Introduction 

The fact that the interstellar medium (ISM) has a hierarchical and self-similar structure 
when observed with sufficiently high dynamic range, is interpreted as evidence of an under¬ 
lying fractal structure (Scalo 1990). The boundaries of the projected images of interstellar 
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clouds are irregular curves whose fractal dimension is around ~ 1.3 (e.g. Falgarone et al. 
1991; Lee 2004). Apparently this is a universal result which does not depend on whether 
tracers of atomic, molecular or dust components are used, whether clouds are selfgravitating 
or not, etc (Williams et al. 2000). The fractal dimension of the projected boundaries is usu¬ 
ally associated with a three-dimensional fractal dimension Df ~ 2.3 for the ISM (e.g. Beech 
1992). It has been argued that a fractal ISM with Df ~ 2.3 could account for the observed 
mass and size distributions of the interstellar clouds (Elmegreen & Falgarone 1996), as well 
as for the intercloud medium properties (Elmegreen 1997a), and even for the stellar initial 
mass function (Elmegreen 1997b, 1999). 

In a previous paper (Sanchez et al. 2005) we studied the effect that the projection of 
clouds has on the estimation of the fractal dimension, and we concluded that a value around 
~ 1.3 for the projected boundaries is more consistent with three-dimensional clouds having 
Df ~ 2.6±0.1. The application of A-variance techniques to Polaris Flare cloud by Stutzki et 
al. (1998) yielded a fractal dimension for the cloud surfaces ~ 2.6 (see also Bensch et al. 2001). 
Elmegreen (2002) simulated fractal brownian motion clouds with average fractal dimension 
~ 2.75 and obtained properties in gross agreement with observations. A fractal medium 
with a relatively high Df value can reproduce observations of HII regions surrounding stars 
(Wood et al. 2005). Henriksen (1986, 1991) used a gravitationally driven turbulence model 
to study the properties of giant molecular clouds suggesting that a dimension ~ 2.7 could 
be necessary to explain the observed properties. Also Fleck (1996) analyzed the properties 
of the turbulent, non-self-gravitating, neutral component of the ISM by using a model of 
compressible turbulence, concluding that the compression parameter that better reproduces 
observations is such that Df ~ 2.5. 

It is important to quantify the degree of complexity (through, for example, the fractal 
dimension) as a first step towards understanding the physical mechanisms responsible for 
structuring the ISM. In models of fractally homogeneous turbulence the fractal dimension of 
fully developed turbulence is 2.5 < Df < 2.75 (Hentschel & Procaccia 1982). On the other 
hand, turbulent diffusion in a incompressible medium generates structures with Df ~ 2.3 
for a Kolmogorov spectrum (Meneveau & Sreenivasan 1990). It is not obvious, however, 
that the energy can cascade down without any dissipation or injection, and deviation from a 
Kolmogorov spectrum is, in general, expected in the ISM (Brunt & Heyer 2002a,b). The ISM 
is a highly compressible and turbulent medium and its fractal dimension depends, among 
other factors, on the degree of compressibility (Fleck 1996) or on the Mach number (Padoan 
et al. 2004). Also, self-gravity could by itself explain many observed properties (de Vega et 
al. 1996), although this fact does not imply that turbulence is not an important factor in 
the ISM. 
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Here we are interested in understanding the relationship between the physical properties 
of the interstellar clouds and their fractal structure, and in verifying whether observed prop¬ 
erties are in agreement with relatively high fractal dimension values (as suggested in Sanchez 
et ah 2005). In other words, we investigate what physical properties are intimately connected 
to the cloud “geometry” and whether this geometry, no matter how it was originated, de¬ 
fines its next evolutionary stage. Our approach in this work is to calculate and analyze the 
properties resulting from a hierarchical structure with a known and perfectly defined fractal 
dimension, no mattering the physical processes behind this structure. In order to do this, 
we simulate interstellar clouds by using a simple algorithm explained in section 2, which 
generates a distribution of points with a very well defined fractal dimension. In section 3 we 
address the density fields resulting from the generated fractals and in section 4 we discuss 
the different ways to construct “clumps” from these density fields. Section 5 is devoted to 
study the properties (masses, densities, sizes, and mass and size distributions) of the derived 
clumps and to compare these results with observations. Finally, the main conclusions are 
summarized in section 6. 


2. Simulated fractal clouds 

To generate fractal clouds we have used the same procedure described in Sanchez et al. 
(2005). Within a sphere of radius Rf we randomly place the centers of N spheres of radius 
Rf/L with L > 1. In each sphere we again place the centers of N smaller spheres with radius 
Rf/L 2 , and so on, up to a given level H of hierarchy. At the end of this procedure there are 
N h points distributed in the space with a fractal dimension given by Df = log(iV)/ log(L). 
This kind of fractals mimics in some way the hierarchical fragmentation process occurring 
in molecular cloud complexes. We have used N = 3 fragments through H = 9 levels of 
hierarchy with the fractal dimension ranging from 2.0 to 3.0. To reduce possible random 
variations all the properties shown are the result of calculating the average of 10 different 
realizations (random fractals). In order to prevent the appearance of a multifractal behavior 
it is necessary that spheres do not superpose when generating the fractal cloud, and this 
requirement prevents the algorithm from generating random fractals with Df > 2.6. Fractals 
with Df = 3 have been obtained by distributing randomly N H points in the available volume, 
and in this case we have calculated the average properties of 50 different realizations. 

The properties resulting from doing an ideal random sampling throughout the hierarchy 
in these kind of fractals have been discussed by Elmegreen (1997a). Let us call Rf and Mf 
the radius and mass (number of points) of the fractal, respectively. The average density 
of the whole fractal is pf = Mf /(4/3)irR 3 f. In the i-th level (i = 0,1,...,77) there are 
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n t = N‘ l fragments with radii RJRf = N~ l ^ D R masses Mi/Mf = N~ l , and average densities 
Pi/pf = jVd 3 /- 0 /^ 1 ). Then, the mass increases with the radius as Mi/Mf = ( Ri/Rf) D f and 
the density varies as Pi/pf = (Ri/Rf)~^~ Df \ The number of fragments with a given mass 
value is n(M) = (M/Mj) -1 , so that the mass distribution function in equal intervals of M 
is dn(M)/dM ~ M -2 , independent of the fractal dimension Df. This means that the mass 
distribution function per logarithmic mass interval is dn(M)/d logM ~ M _1 . On the other 
hand, the number of fragments depends on the radius as n(R) = (R/Rf)~ D f and the size 
distribution function per logarithmic interval will be dn{R)/d\ogR ~ R~ D f. Therefore, the 
theoretically predicted mass and size distributions from random sampling in a fractal cloud 
are power laws of the form dn(M)/d\ogM ~ M~ aM and dn(R)/d\ogR ~ R~ aR with indices 
(i.e., slopes in a log-log plot) given by = 1 and an = Df , respectively. However, in the 
hierarchical structure the smaller fragments are not isolated objects but they are included 
into the big ones, and if double counting is avoided then the indices and a r become lesser 
(Elmegreen 1997a). Even though we do not take this effect into account, smaller fragments 
can be blended into big ones in real clouds. This blending can be “real” (random cloud-cloud 
coalescence) or “artificial” (due to resolution limitations), but the effect should be the same: 
to decrease the expected values of oim and or. Therefore, it is not clear what values should 
we measure when observing fractal complexes as the ones described above, and this is the 
main goal of this work. 


3. Probability density function 

Before proceeding with the calculation of the cloud properties, we have to estimate 
the one-point probability density function (hereinafter pdf) of the density for the fractal 
distribution of points. This is an important issue because we will need the density pdf to 
compute clump masses and radii in the following sections. The most direct way to do this 
is to place a grid of cells on the fractal, then calculate the density in each cell and plot a 
histogram with the number of cells in a certain range of density. An alternative way that 
guarantees a good sampling is to place randomly a high enough number of cells within the 
fractal. The density in each sampling is then the number of points inside the cell divided 
by the volume. As critical requirement for a suitable estimation of the pdf, the cell size 
(. R c ) has to be big enough to produce a large number of density values, i.e., it has to be 
much bigger than the minimal distance between two points in the fractal. Figure 1 shows 
the density pdf (per logarithmic interval of density) resulting from using spherical cells with 
radius R c = -R//10 for the fractals with dimensions Df = 3 and Df = 2.6. For the case 
Df — 3 (very homogeneous distribution of points) we see that, as expected, the density pdf 
is very similar to a gaussian with the maximum around p/pf = 1. For comparison we also 
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p/pf 

Fig. 1.— The density pdf calculated in a direct way by using cells with radius R c = Rf /10 
for the fractals with dimensions Df — 3 (dashed line histogram) and Df — 2.6 (solid line 
histogram). The dashed curve is a gaussian fit for Df = 3 and the solid curve shows the 
density pdf calculated for Df = 2.6 but using a gaussian kernel with a = R c (see text). 
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show, with a dashed curve, a gaussian distribution fitted to the Df — 3 histogram. For 
the case Df = 2.6 the resulting pdf departs from a gaussian distribution. The behavior 
resembles a power law at low densities (with slope ~ 0.65) with a sharp cnt-off after a 
maximum. The expected density value, when a whole structure with size R c is sampled, is 
PciPf = {Rc/ RfY^~ Df \ which for Df = 2.6 means pj Pf — 2.5. The pdf cut-off is close 
to this value, although higher values can be seen due to the presence of random denser 
regions. However, much higher density values cannot be detected by using this cell size, it is 
necessary to use smaller cells to be able to sample efficiently denser structures belonging to 
lower levels. The average density of the (upper) level where the cell of size R c is embedded 
is given by p c _i/p/ = (N l ! D fRJRf)~^~ D f\ which gives ~ 2.12 for Df = 2.6. This is indeed 
the value around which the maximum of the distribution is found. Smaller density values 
can be obtained as the (fixed volume) cell samples less massive structures in lower levels. It 
has to be mentioned that cells of different sizes tend to produce distributions with different 
positions for the maximum (higher values as R c decreases), but the distributions always 
follow roughly a power law at low densities with a cut-off after the maximum. 

The pdf estimated in this direct way exhibits a discrete distribution because there can 
only be an entire number of points inside the cell. In order to turn it into a continuous 
distribution we use a kernel function to convolve the data and calculate the density through 
the whole available volume. We have adopted a gaussian kernel (normalized to have an 
integral of 1); thus, the density at the position ( Xi,yi,Zi ) is (Silvermann 1986) 


/ i x i, Vi, Zi) 


1 

N H a 3 {2ix ) 3 / 2 


n h 

Y (, xp 

3 = 1 



- Xj ) 2 + {yi - Vj ) 2 + {Zi 
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where o is the “smoothing” parameter, which determines the volume used to estimate the 
average density in each position (roughly equivalent to the former cell size). As an example, 
Figure 1 shows with solid line the pdf calculated for Df = 2.6 by using equation (1) with a = 
R c . As before, the pdf was constructed by calculating the local density in randomly chosen 
positions and plotting the number of density values per logarithmic interval of density. The 
result is very similar to that from the direct method except that now we have a continuous 
distribution. Again, the exact functional form will depend on the adopted values for cr. 
Large a values (of the order of the fractal size) lead to a too homogeneous distribution of 
density values, whereas small values (much smaller than the minimal distance between two 
points) tend to produce very narrow density peaks. To avoid both undesired situations and 
to be as objective as possible we always have used the cr values that maximize the likelihood 
(Silvermann 1986), that we call a opt herein. Figure 2 is an example of the appearance of the 
density field resulting from applying equation (1) to a fractal with Df = 2.6. Each map is 
a slice along the z = 0 plane showing density contours of a part of the three-dimensional 
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Fig. 2.— A slice along the z = 0 plane showing density contours of part of a fractal with 
Df = 2.6, calculated by using equation (1) with (a) a = a opt and (b) a = 2 a op t- 



density field. The left-side map (a) was generated with a = a opt , whereas the right-side 
map (b) is exactly the same region but generated with a = 2a opt . Clearly we see that as 
a increases some density peaks tend to blend with other peaks forming smoother regions. 
Several tests showed that the properties studied in this work do not depend (within the error 
bars) on the exact value of a as long as this value is of the order of a op t • This value can be 
different from fractal to fractal but always was of order 10 ~ 2 R f . 

In Figure 3 we show the normalized density pdf calculated (with a = a opt ) for the 
fractals with dimensions Df — 3, Df — 2.6, Df = 2.3, and Df = 2.0. The long tails at low 
densities come from the sampling in or close to the edge and/or the voids existing in the 
fractal clouds. The maximum of the distribution moves toward higher density values as Df 
decreases according to the relation p max oc R-( 3 ~ D f'> for R ~ constant , from p max /Pf — 1 
for Df = 3 to Pmax/Pf — 10 for Df = 2.0. There is always a turn over in the pdf, decaying 
rapidly at densities just after the maximum. 

The density pdf is a fundamental statistical quantity characterizing the structure of a 
given medium. Despite the pdf being a one-point statistic containing no spatial information, 
the knowledge of this function is a fundamental step prior to the analysis of the star formation 
process and the stellar initial mass function (e.g. Padoan & Nordlund 2002; Scalo et al. 
1998). The problem is that the density pdf cannot be directly measured in the ISM, instead 
the column density pdf is measured. Ostriker et al. (2001) used three-dimensional numerical 
simulations of isothermal magnetohydrodynamic turbulence to compare the resulting density 
field with the column densities, hireling that both distributions had very similar shapes 
(approximately lognormal). Also Fischera & Dopita (2004) showed that lognormal density 
distributions produce lognormal column density pdf’s. However, column density observations 
cannot always be fitted to lognormal distributions and can show extended tails (Blitz & 
Williams 1997) or simply a power law behavior (Heilcs & Troland 2005). 

The suggestion that the density pdf in the ISM should be consistent with a lognormal 
distribution comes from the central limit theorem applied to a multiplicative hierarchical 
density held (Vazquez-Semadeni 1994). Several authors have found that the density pdf of 
turbulent gas should exhibit a roughly lognormal distribution under a wide range of physical 
conditions (Ballesteros-Paredes 2004; Elmegreen & Scalo 2004, and references there in). For 
a polytropic gas (where P oc p n ) the shape of the pdf depends on the effective polytropic 
index n. If the gas is isothermal (n = 1) the pdf is lognormal (Vazquez-Semadeni 1994) but 
extended tails can develop at high or low densities if n ^ 1 (Passot & Vazquez-Semadeni 
1998; Scalo et al. 1998; Ostriker et al. 2001), and the departure from a lognormal is larger 
at higher Mach numbers (Beresnyak et al. 2005). The one-dimensional simulations of highly 
compressible turbulence of Passot & Vazquez-Semadeni (1998) showed that a power law tail 
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P/Pf 


Fig. 3.— The density pdf calculated by using a gaussian kernel with o = a opt , for different 
values of fractal dimension: Df = 3 (dotted line), Df = 2.6 (solid line), Df = 2.3 (long 
dashed line), and Df = 2.0 (short dashed line). 
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develops at relatively high densities if n < 1 and at low densities if n > 1. Their results for 
n > 1 and high Mach numbers resemble the pdf’s we showed in Figure 3. 

The situation becomes more complicated if we try to understand the dependence of the 
density pdf on the physical mechanisms acting in the turbulent medium. Li et al. (2003) 
found that shocks sweep up the gas producing low density regions and pdf’s skewed to low 
densities, but self-gravity tends to collect the gas producing pdf’s with positive skewness. 
The presence of a magnetic field has, however, very little effect on the density pdf (Passot 
& Vazquez-Semadeni 2003). The large-scale simulations performed by Wada & Norman 
(2001) showed a roughly constant pdf shape (lognormal at high densities and normal at low 
densities) regardless of whether star formation is considered or not in the simulation. Slyz 
et al. (2005) arrive to very different conclusions, finding that the pdf does depend sensitively 
on the simulated physics. Their pdf is consistent with a lognormal function only for the 
runs where neither supernova feedback nor self-gravity are implemented. With the addition 
of self-gravity the pdf approaches a power law at high densities, but it becomes markedly 
bimodal when supernova feedback is included (with or without self-gravity). 


4. Clump definition 

Starting from the three-dimensional density held calculated by applying equation (1) 
to the generated fractals, we want to study the properties of the resulting “density con¬ 
densations”. We will call these regions “clumps” to differentiate it from the whole fractal 
structure, which we could associate with a giant molecular cloud. Actually this denomination 
is arbitrary, firstly because there are no precise definitions of terms like “cloud”, “clump”, 
or “core” (Larson 2003), and secondly because here we can not speak about an absolute 
spatial scale but about self-similar structures inside a region of size Rf. Different schemes 
can be adopted in defining a clump. For example, the algorithm called GAUSSCLUMPS 
(Stutzki & Gusten 1990) uses a least squares fitting procedure to decompose iteratively the 
emission in one or more gaussian clumps, instead CLLIMPFIND (Williams et al. 1994) asso¬ 
ciates each local emission peak and the neighboring pixels with only one clump (similar to 
the usual eye-inspection procedure). In both cases, the implicit assumption that the radial 
velocity coordinate can be replaced by the radial distance is made, but this assumption is 
not necessarily always satisfied (Ostriker et al. 2001). In our case, this problem does not 
exist because we are using three spatial coordinates and density values, but we have to keep 
in mind that comparison with observations may not be done directly. For simplicity, we have 
chosen to associate each peak in density with one clump (in the same way that CLUMPFIND 
does). Once we have identified a peak we determine the mass of the clump by integrating 
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the density field over the neighboring region. In order to proceed with the calculation we 
need to place a grid over the volume occupied by the fractal cloud. The size of the cell was 
chosen to equal de value of a opt , that is, to equal the resolution with which the density field 
was generated. Since a opt /Rf was always of order of 1CD 2 or less, there were always ~ 10 6 
cells or more over the fractal. 

An important aspect to consider when defining a clump is the noise present in the data, 
whether it be observational or numerically simulated. Usually the problem is approached 
by defining a free parameter such that fluctuations of the order of (or smaller than) this 
parameter are not considered as clumps. For example, in the algorithm CLUMPFIND 
this parameter is the step between two successive contours, in the simulations of Gannnic 
et al. (2003) is their smoothing length, and in this work it is the parameter a used in 
equation (1) to generate the density field. Another very important point is to define the 
boundary of the cloud, i.e., to define which neighboring region “belongs” to the clump and 
which not. Again the usual approach is to introduce an additional parameter: the threshold 
density (or emission). Thus, a common definition for clump in literature is all the pixels 
around a maximum in density with values greater than some predefined threshold value. 
The algorithm we have constructed to identify clumps works in the following manner. First, 
it associates the pixel with highest density value to the first clump. Then it takes the pixel 
whose density value is just below the highest, and if this pixel is neighbor of the first one 
then it “belongs” to the first clump. If by the contrary this pixel is isolated then it is labeled 
as the first pixel of the second clump. All the pixels are explored successively in decreasing 
order of density in the same way: if a pixel is contiguous to the boundary of some existing 
clump then it belongs to this clump, otherwise it belongs to a new clump. This procedure 
continues until some threshold density is reached. 

When this threshold is relatively low there can be an additional problem: two (or more) 
peaks in density may be present in the region above the threshold, and we have to “decide” 
to which of the available peaks the pixels that are contiguous to two or more boundaries 
should be assigned to. The criterion normally used is purely geometric: each “ambiguous” 
pixel is assigned to the clump with the closest peak (just like CLUMPFIND does). However, 
we should ask ourselves whether different criteria could be used to assign “membership” in 
a more “suitable” way. An alternative criterion could be to assign each ambiguous pixel to 
the clump with highest “gravitational potential”. We have tested the former criterion by 
estimating the potential (at the moment of considering the i-th pixel) as the current clump 
mass (Mj) divided by the peak-pixel distance (d*). Then, each ambiguous pixel is attached 
to the clump with higher Mj/dj value rather than to the clump with higher 1/dj value (i.e., 
smaller distance). If the threshold density is close to the maximum density values for two 
neighbor clumps, then the clumps will be physically separated and any criteria for ambiguous 
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pixels will yield the same result. As threshold density (pth) decreases the adjacent surface 
between two neighbor clumps increases and therefore the effect of different criteria becomes 
more important. As an example, Figure 4 shows ten of the resulting clumps when p t ^ = 0 for 
the first random fractal we simulated with Df — 2.6 and for both criteria mentioned above. 
The criterion of distance tends to produce more symmetrical (and similar in shape) clumps, 
and in this case the clump volume distribution relates directly with the peak-peak distance 
distribution between neighbors, which is random for the fractal sets we are generating. This 
degenerates into a clump mass distribution narrower than for the criterion of potential, as can 
be seen in Figure 5a where the number of clumps per logarithmic mass interval (normalized 
to the total number of clumps) N ci /N tot has been plotted as a function of the clump mass in 
units of the total mass ( M c i/Mf ). Thus, the same cloud density held and therefore the same 
number of clumps yields very different mass distributions only changing the criterion with 
which the low density pixels are distributed. The differences are amplified simply because we 
are artificially creating clumps from extremely low density material which typically remains 
in the form of inter-clump gas. The clumps obtained by distributing all the cloud mass 
in clumps are far from looking similar to those obtained from observations: there are no 
observed clouds formed only by clumps without inter-clump material. We have verified that 
in more realistic situations when p t h > 0, i.e. when an inter-clump density is defined, clumps 
resulting from applying different criteria and therefore their global properties (such as mass 
and size distributions) gradually begin to converge as p t h increases. In fact, when only 10% 
of the cloud mass is forming clumps both criteria for selection of clumps yield very similar 
mass distributions (see Figure 5b). In other words, when we only see the dense “cores” the 
method of distributing low density pixels does not play any important role. In section 5.3 
we analyze the effect of the threshold density on the global clump properties in order to 
understand how these properties are modified in clouds observed with different degrees of 
sensitivity and contrast. 

We have tested several different criteria for clump selection but the global results always 
are similar to one or another of the two criteria above, depending on whether peak-pixel 
distance or clump mass is the dominant factor in the criterion applied. Here we have chosen 
to discuss in detail the properties resulting when using the (physically inspired) criterion 
of the “potential”. When considering the clump mass for deciding about ambiguous pixels 
this criterion takes into account in some way the underlying density pdf, which depends on 
the cloud fractal dimension (section 3). The properties we derived by using the criterion of 
distance map the random distribution of distances between neighbor density peaks rather 
than of intrinsic density structure. However, the sensitivity of the cloud properties, both 
simulated and observed, to different clump selection criteria is an important problem that 
must be addressed in future works. 
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Fig. 4.— Images of the clumps resulting (each one with a different color) from ten density 
peaks when p t h = 0, for the first random fractal generated with Df = 2.6, using (a) the 
criterion “potential” and (b) the criterion “distance” (see text). The coordinates are the 
pixel coordinates. 



Fig. 5.— The clump mass spectra for ten random fractals with Df = 2.6 when using 
the criterion of potential (solid line histograms) and the criterion of distance (dashed line 
histograms), (a) All the cloud gas is forming clumps and (b) only 10% of the cloud gas is in 
the form of clumps. 
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5. Properties of fractal clouds 

We have defined each clump as all the mass enclosed from the density peak to a threshold 
density. For convenience, we have chosen the simplest way to define this threshold: a constant 
value such that a given ratio, e, of the total mass in clumps to the total cloud mass is 
obtained. If e = 1 the number of clumps formed will equal the number of density peaks, 
which is given by the three-dimensional density field (i.e., by Df) and it does not depend on 
the clump selection criterion. Low Df values imply more fragmented clouds (higher number 
of clumps). But additionally at low e values some peaks could be below the threshold density 
and then there could be a lesser number of clumps. Table 1 shows the average number of 
clumps (per simulation) obtained for the results we discuss in this work. We first discuss the 
results for e = 1 and in section 5.3 we analyze the effect of changing the threshold density. 


5.1. Masses, densities, and sizes 

Generally the resulting clump shapes do not have a regular geometry (see the examples 
in Figures 2 and 4). Linder these circumstances it is not easy to define a clump “radius”, 
so we have defined a characteristic radius simply as the cubic root of the clump volume. 
Figure 6 shows the mass of the clumps ( M c i ) as a function of the radius ( R c i ) for four 
different values of the fractal dimension (Df = 2.0, 2.3, 2.6, and 3.0). The corresponding 
average densities (p c i) are shown in Figure 7. The range of possible sizes depends in principle 
on the fractal dimension (that determines the way in which fragmentation occurs) and on the 
number of levels (that determines how many fragmentations occur). When generating the 
fractal structure the smallest fragments have radii R m i n /Rf — (1 /N H ) l ^ D f. Computational 
limitations prevented us from using more than H = 9 levels, so the smallest structures have 
sizes ~ 7 x 10 ~ 3 Rf for Df = 2.0 and ~4x 10 ~ 2 i?/ for Df = 3.0. The size range obtained 
is more or less between these minimum values and the theoretical maximum R c i/Rf = 1. 
Generally speaking, mass does not obey a power law with size of the type M c i oc RR d along all 
the size range (Figure 6 ), and therefore neither density obeys the corresponding p c i <x i ?^ -3 
(Figure 7). We would expect 7 = Df in the case of random sampling throughout the 
hierarchy (section 2), but what we see is 7 ~ 3 (i.e., p ~ constant ) at small sizes and 
then 7 gradually decreases as R c i increases. For comparison, we show in Figure 6 the lines 
corresponding to p c i = constant. The decrease in the 7 value is more abrupt for low Df 
values, while for the extreme case Df = 3 we get 7 ~ 3 along almost the full size range. For 
Df = 3 the range in density values is very narrow and close to p c i/Pf = 1 (as expected), 
but for lower Df values the density range increases and tends to higher values (Figure 7 ). 
The latter situation corresponds to more fragmented structures with small dense clumps 
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Table 1. Average number of clumps formed per fractal cloud 


Df 

e = 1 

e = 0.5 

e = 0.25 

e = 0.1 

2.0 

764.6 

764.6 

737.6 

556.8 

2.1 

620.6 

620.6 

592.7 

469.3 

2.2 

531.1 

531.1 

514.6 

410.2 

2.3 

427.5 

427.2 

412.5 

316.8 

2.4 

343.9 

343.4 

328.7 

244.2 

2.5 

262.1 

262.0 

243.5 

175.4 

2.6 

189.8 

189.5 

175.0 

133.6 

3.0 

118.0 

116.4 

106.1 

80.4 



R c! / Rf R c! / R f 


Fig. 6.— Clump mass, M c [, as a function of radius, R c i, for the labeled values of fractal 
dimension Df. The solid line corresponds to the constant density case. 
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Fig. 7.— The average clump density, p c i , as a function of the radius, R c i, for the labeled 
values of fractal dimension Df. 




separated by large low density regions. 

An empirical power law scaling relation with 7 ~ 2 was first pointed out by Larson 
(1981) for molecular cloud complexes (in the range 0.1 < R < 100 pc and 1 < M < 10 6 M 0 ), 
although it has been argued that this result could be simply an artifact due to observational 
limitations that strongly constraint the range of observed column densities (Scalo 1990). 
This power law behavior has been generally interpreted as a consequence of the mechanical 
equilibrium in self-gravitating, turbulent molecular clouds, but also very different physical 
processes can reproduce this type of relations (see Elmegreen & Scalo 2004, and references 
therein). From the observational point of view, the problem is that many studies show a 
weak or scatter-dominated correlation. Moreover, comparing different observational studies 
is always difficult because of the existing variety in observational techniques and criteria used 
to define or derive properties such as mass, radius, etc. Elmegreen & Falgarone (1996) used 
several surveys from the literature and found 7 ~ 2.35, which was interpreted as a direct 
consequence of an underlying fractal structure with Df ~ 2.35. However, this result only 
emerges when the clouds and clumps from many regions are analyzed as an ensemble, while 
the power laws measured in individual regions have 7 ~ 2.4 — 3.7, with values > 3 being 
typical. Reid & Wilson (2005) found 7 ~ 1.5 — 2.1 in the range 1 < M < 10 3 M 0 , whereas 
Heithausen et al. (1998) found a steeper slope (7 ~ 2.3) in the range 10~ 4 < M < 10 2 M 0 . 
Caselli & Myers (1995) have already observed a change in the power law slope between large 
mass and low mass cores in Orion A and B. Falgarone et al. (2004) have argued that, in 
spite of the large scatter, it seems that 7 ~ 2 better characterizes large scale structures and 
7 ~ 2.3 the small scale ones. 

Our results seem to favor a change in the slope of the M d — R d relation, although 
less steep as Df increases. The manner in which M d increases with R c i is determined by 
the density profile of the clumps. For flat density profiles, p(x, y, z ) = constant, we would 
expect to see M d = f y ^ p(x, y, z)dV ~ R d . For extremely narrow density profiles of the form 
p(x, y, z) = poS(xo, yo, zo), being 5 the delta of Dirac and po the density value of a little region 
of volume Vo around (xo,yo,zo), we would have M d = f v ; p(x,y, z)dV = poVo = constant. 
The observed behavior is directly related to a gradual change in clump density profiles from 
very flat shapes at small sizes to narrower ones at large sizes. As an example, Figure 8 shows 
density profiles for two different clumps in one of the fractals with Df = 2.6. We have plotted 
the density of each pixel as a function of its distance from the pixel having the maximum 
density value. The scatter in these plots is clue to the asymmetry of the clump shapes. The 
figure labeled (a) corresponds to a clump ~ 40 times more massive and ~ 10 times larger 
than the clump in (b), and the characteristics mentioned before can readily be appreciated. 
This is an interesting property of the fractal density distribution we are generating, in which 
large clumps are the result of the “agglomeration” (through a gaussian convolution with 



ity of e; 
ng the 
Is with 









equation 1) of particles in the last level of hierarchy in the original fractal. Its consequences 
for the derived clump properties will be discussed in section 5.3. 


5.2. Clump mass and size spectra 

The clump mass spectra (number of clumps N d per logarithmic mass interval) for differ¬ 
ent Df values are shown in Figure 9. The filled histograms indicate the ranges corresponding 
to clumps having a volume larger than ~ 10 pixels, which were the ranges where a power 
law-like behavior was observed: N c j oc MJ“ M . The best fits in these ranges (solid lines) show 
clearly steeper slopes as the fractal dimension increases. The fits yielded: ajf ~ 0.33 ± 0.12 
for D f = 2.0, a M - 0.36 ± 0.08 for D f = 2.3, a M - 0.54 ± 0.06 for D f = 2.6, and 
cxm — 0.64 ± 0.03 for Df = 3.0. The mass spectra become flatter in the low mass range. 
We have to point out that we are not saying that the intrinsic mass spectra are actually 
power laws, but we are simply characterizing the spectra with the functional form commonly 
used by observers in order to better compare these results. From simple random sampling 
through the hierarchy we would expect = 1 (see section 2), but we already mentioned 
that if double counting is avoided and/or random blending is taken into account the observed 
index should be less than this value. Clouds with high fractal dimensions arc less fragmented 
than low fractal dimension clouds. They have a relatively small number of clumps but, on 
average, these clumps have larger sizes and higher masses. As Df increases clouds tend to 
be more homogeneous, and that is why the average clump densities are always close to the 
whole cloud density (Figure 7). Additionally, the volume in clumps approach the total cloud 
volume as Df increases (the filling factor tends to 1 as Df tends to 3), and the clumps are 
relatively close to each other. On the opposite, at very low fractal dimensions the high den¬ 
sity clumps are separated by low density (or empty at all) regions. Ultimately, low fractal 
dimension clouds distribute their material more homogeneously between the larger number 
of well differentiated clumps. At high fractal dimensions a greater amount of “ambiguous” 
material is assigned to massive clumps leaving a higher number of small clumps. It is for 
this reason that the mass spectrum slope is flatter for small Df values. 

The clump size spectra (per logarithmic size interval) are shown in Figure 10 for the 
same fractal dimension values as in Figure 9. Again, the filled histograms indicate the ranges 
for clumps larger than ~ 10 pixels, and the solid lines are the best fits in these ranges. The 
near-power law behavior N d oc 7?^"^ becomes steeper as Df increases, as the fitted slopes 
show: a R ~ 0.99 ± 0.07 for D f = 2.0, a R ~ 1.15 ± 0.08 for D f = 2.3, a R ~ 1.42 ± 0.13 for 
Df = 2.6, and a R ~ 2.09 ±0.08 for Df = 3.0. As for the mass spectra, the final distribution 
is flatter than the random sampling slope which is a R = Df and, as before, the relative 
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M C | / M f M c , / M f 

Fig. 9.— The clump mass spectra for the labeled values of fractal dimension Df. The filled 
histograms indicate the mass range for clumps having a volume larger than ~ 10 pixels, and 
the solid lines are the best fits in these ranges. 
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Fig. 10.— The clump size spectra for the labeled values of fractal dimension Df. The filled 
histograms indicate the size range for clumps having a volume larger than ~ 10 pixels, and 
the solid lines are the best fits in these ranges. 
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importance of the clump selection criterion when applied to more or less fragmented clouds 
causes the differences in the derived slopes for different Df values, ft is interesting to mention 
that if mass and radius could be related as M d oc F0 d then the power law indices satisfy 
a relationship of the form an — 7 olm- For the random sampling case we have 7 = Df, 
otM = 1, and then an = Df, as shown in section 2; but our results do not show a scale 
relation between mass and radius (7 ^ constant, see Figure 6) and therefore this kind of 
simple relations cannot be found. 

The power-law behavior obtained for both the mass and size distributions are in gen¬ 
eral consistent with observations. However, the detailed comparison is extremely difficult 
considering the great diversity of results reported in the literature and considering that no 
clear relationship has been observed between the mass or size distributions and the mean 
physical properties of the ISM. The most accepted view is a nearly constant mass and size 
distributions throughout a large range of scales, according to a fractal picture for the ISM 
(see, for instance, the reviews of Evans 1999; Williams et al. 2000; Falgarone et al. 2004, 
and references therein). When a sufficiently large number of clumps is considered, galactic 
surveys indicate that on average q;m ~ 0.8 in the range 1 < M < 10' M 0 and an ~ 2.3 in 
the range 0.1 < R < 10 2 pc (Elmegreen & Falgarone 1996), and these results seem to depend 
very little, if any, on physical properties like selfgravitation or the degree of star formation 
activity (Simon et ah 2001). The extension at lower masses and sizes yields 0.6 < 0.8 

for 10~ 4 < M < 10 3 M q and an ~ 2.0 ± 0.2 for 10~ 2 < R < 1 pc (Heithausen et ah 
1998; Kramer et ah 1998). The mass and size distributions obtained here have indices in 
the range 0.3 < < 0.6 and 1.0 < an 2.1 for 2 < Df < 3. If we assume Df ~ 2.6 

then a,M ~ 0.5 and an ~ 1.4, a little lower but in gross agreement with mean observed 
values. Steeper slopes can be obtained only for Df > 2.6. These results seem to favor the 
idea that the ISM has a fractal dimension higher than the usually assumed value Df ~ 2.3 
(Sanchez et ah 2005). However, the situation is far from being totally understood. The 
mass distribution slopes may be as steep as, for example, ay,/ ~ 1.1 for filamentary clumps 
in Orion A molecular cloud (Nagahama et ah 1998) or as flat as qm — 0.3 for the Rosette 
molecular cloud (Williams et ah 1994), not being obvious whether these observations reveal 
real differences from region to region. Schneider & Brooks (2004) have demonstrated that 
the derived properties differ significantly when using different clump identification meth¬ 
ods. The algorithm GAUSSChUMPS (Stutzki & Gusten 1990) is able to separate blended 
clumps, although tends to introduce spurious features if the structures do not have gaussian 
shapes. On the other hand, CLUMPFIND (Williams et ah 1994) is able to find arbitrarily 
shaped clumps but it does not detect less massive clumps which are blended with others. 
Thus, generally speaking, one would expect steeper mass distributions for GAUSSCLUMPS 
than for CLUMPFIND (Schneider & Brooks 2004), but error bars are usually large and this 
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behavior cannot be appreciated (Mookerjea et al. 2004). The situation is worst for the size 
distributions because clump sizes are quantities more difficult to estimate in a robust man¬ 
ner, particularly for small, low intensity objects close to the resolution limits (Elmegreen & 
Falgarone 1996; Heithausen et al. 1998). 


5.3. Dependence on the threshold density 

The previous results refer to the case e = 1. A smaller e value implies that not all 
the cloud complex mass has formed clumps. We have simulated this effect by considering 
a (constant) threshold density such that pixels with densities below this threshold value 
are not taken into account when constructing clumps. In principle, e could be related to 
observational limitations, because smaller e values arise from the fact that there is a threshold 
density level for observations. Also we could associate it with the actual mass fraction in 
clumps, being the rest of the mass kept in a diffuse interclump medium (Blitz & Stark 1986; 
Wood et al. 2005). In such a case e could be around ~ 0.9 or even much less (Blitz & 
Stark 1986; Pagani 1998). One possibility is to assume that a fixed fraction of clump masses 
will form dense prestellar cores and then e ~ 0.1 (if we assume a star formation efficiency 
~ 10%), although obviously the formation of gravitationally bounded cores is a process more 
complex than this. The threshold densities we discuss here have been chosen to obtain the 
following e values: 0.5, 0.25, and 0.1, and all the clump properties were recalculated. 

For e < 0.5 the mass-radius relations for the clumps correspond to the p c i = constant 
case. The range of masses decreases, which is not a surprising result since part of the clumps 
mass has been removed, but additionally the mass spectra behavior modifies in an interesting 
way: as e decreases the slope in the power law range becomes steeper. The mass spectra 
for the same fractal dimensions as before but for e = 0.1 are shown in Figure 11. The same 
general behavior is kept (a power law at high masses with a flattening at low masses) but the 
slope values changed notoriously. For example, for the case Df = 2.6 the index increased 
from 0.54 ± 0.06 (e = 1) to 1.19 ± 0.20 (e = 0.1). We have summarized all the obtained 
mass spectra by plotting in Figure 12 the power law slope cxm as a function of the fractal 
dimension Df for three different values of e (1.0, 0.25, and 0.1). The result for e = 0.5 is so 
close to the result for e = 1 that it is not shown for simplicity. As already mentioned, the 
general behavior is an increase of cc M (steeper slopes) as Df increases, but also increases 
as e decreases. In spite of the bar sizes, Figure 12 allows us to quantify approximately the 
dependence of q?m on both Df and e. Part of the uncertainty comes from fitting power law 
functions to distributions that depart from this behavior. 

Concerning the clump size spectra, when e = 0.5 the size range decreases in such a way 





Fig. 12.— The power law index chm as a function of the fractal dimension Df for three 
different values of e: 1.0 (circles), 0.25 (squares), and 0.1 (triangles). The bars on the data 
are the standard deviations of the best fits. The results for e = 0.5 are very close to those 
for e = 1 and are not shown for the sake of clarity. 








that a clear power law behavior is not seen, mainly at low fractal dimensions. However, we 
have followed the objective procedure of fitting a power law for structures bigger than ~ 10 
pixels, and we obtained flatter distributions. For example, for Dj = 2.6 we have that olr 
decreased from 1.42±0.13 (e = 1) to 0.92±0.20 (e = 0.5). The obtained power law slopes 
as a function of the fractal dimensions Df for e = 1 and e = 0.5 are shown in Figure 13. For 
smaller e values the size ranges were so narrow that fits could not be done and/or standard 
deviation bars were extremely large. The power law index olr increases as Df increases, but 
now we have that a r also decreases when e decreases, in opposite to the mass spectra. 

The trend of steeper mass distributions for higher threshold densities (smaller e values) 
was seen by Elmegreen (2002) in his model of fractal brownian motion clouds, although he 
only simulated clouds with constant fractal dimension (Df ~ 2.75 on average). It is not easy 
to understand the dependence of and on e, but as in the work of Elmegreen (2002) 
it has to do with changes in the density profiles between the smallest and largest clumps. 
In section 5.1 we argued that density profiles of small clumps are flatter than those in large, 
massive clumps (which have lesser average density). Large clumps have narrow, tall profiles 
with long tails at low densities (see the example in Figure 8). Then, when the threshold 
density is increased, both sizes and masses decrease, but in different ways. Large clumps 
decrease their sizes in a higher proportion than the small ones, because of the long tails. 
Thus, the size distribution flattens because the proportion of high radius clumps increases 
(relative to the small ones). On the opposite, there is not too much mass below the low 
density tails and then the masses of large clumps decrease in a lower proportion than for 
smaller clumps. Therefore, when the threshold increases the proportion of high mass clumps 
decreases and the mass distribution becomes steeper. 

We have seen that for more realistic e values (i.e., e < 1) we get steeper mass functions 
and flatter size functions. The case e = 0.5 and Df = 2.6 yields 0.7 < < 1.1, which is 

much flatter than the average observed in molecular clouds but remarkably similar to the 
size distribution recently observed for HI clouds in the LMC (Kim et ah 2005). Additionally, 
according to our results, molecular cloud complexes with Df = 2.6 will have mass functions 
a,M ~ 1.2 ± 0.2 when e = 0.1, very close to the value 1.35 of Salpeter (1955). In fact, the 
observed mass spectrum of dense cores seems to have a form significantly different from 
molecular clumps, with the slope relatively flatter (c— 0.5) at low masses and steeper 
(cxm — 1.5) at high masses (Motte et ah 1998; Reid & Wilson 2005). Our results reproduce 
this behavior, we clearly appreciate (in spite of the large standard deviation bars) a sys¬ 
tematically steeper mass function at high masses as e decreases. This is a very conspicuous 
result: the denser regions in the fractal cloud complex have a Salpeter-like mass function 
steeper than the mass function of the whole complex. An interesting point is that the mass 
distribution slope for e = 0.1 is near to the Salpeter value (within the error bars) in a wide 






range of fractal dimension values (2.0 < Df < 2.6). It seems that the information about the 
fractal density structure is lost when we only see the central dense “cores”. The results we 
showed in this section suggest a direct link between cloud structure and star formation. 


6. CONCLUSIONS 

In this work we have studied the relation between the physical properties of interstel¬ 
lar cloud complexes and their underlying fractal structure. For a better understanding of 
this relationship we have used a simple algorithm that generates fractal clouds with fractal 
dimensions well defined in a wide range of space scales. We have not considered physical 
processes to generate either the fractal structure or the embedded clumps. In any case, the 
exact physical nature of clumps still remains uncertain; whether they are temporary density 
fluctuations caused by supersonic turbulence or more stable structures confined by the inter¬ 
clump medium (Williams et al. 2000). The simple approach given here allows us to analyze 
in an empirical way the dependence of the ISM properties on both the fractal dimension and 
the threshold density. 

We observe that the number of clumps (as given by the number of relative maxima of 
density) depends only on the density structure of the whole cloud, which can be associated 
to a single parameter: the fractal dimension. However, when the whole mass of the cloud 
is distributed amongst the different clumps (i.e., e = 1 ), the mass and size distributions are 
highly dependent on the clump defining criterium chosen. On the other hand, as the fraction 
of the total mass in the form of clumps (e) decreases the mass and size distributions become 
similar for different criteria of clump selection, i.e., when only the “cores” of the “clumps” 
are selected their distributions in mass and radius are not dependent on the selection criteria. 
This fact leads to an interesting conclusion: the “cores” mass and size distributions are only 
driven by the fractal dimension of the clouds. A true molecular cloud contains relatively 
empty voids, dense cores, as well as rarefied intercloud material (Gammie et al. 2003), which 
looks like more similar to the distributions of cores obtained for e = 0.1 than the clumps 
distribution for e = 1 . 

In general, the masses and radii of the resulting clumps do not fulfill simple power law 
relations of the type M d oc fV d along all the size range. We obtain that 7 ~ 3 at small 
sizes and 7 < 3 at larger sizes with the exact behavior depending on the fractal dimension. 
The reason for the departure from a power law has to do with the differences in the density 
profiles between the smallest and largest clumps: small clumps have flat density profiles 
whereas large clumps have narrow and tall profiles with long tails at low densities. 
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The number of clumps per logarithmic interval of mass (or size) can be fitted to a power 
law function A% oc M~ t aM (or A% oc in a certain range of mass (or size). The indices 

% and an depend both on the fractal dimension (Df) and on the fraction of the total mass 
in the form of clumps (e). Our main results are summarized in Figures 12 and 13, and they 
refer to the dependence of the clump mass and size spectra on both Df and e. For the 
case e = 1 we obtain that as Df increases from 2 to 3 aj( increases from ~ 0.3 to ~ 0.6 
whereas increases from ~ 1.0 to ~ 2.1. Rough comparison with observations suggests that 
Df ~ 2.6 is consistent with the average properties of the ISM. This value for the ISM fractal 
dimension is in agreement with previous results based on the fractal dimension measured on 
the projected images of clouds (Sanchez et al. 2005). On the other hand, as e decreases a M 
increases and cxr decreases. For the case e = 0.1 (only 10% of the complex mass is in the 
form of dense clumps) we obtain aM — 1-2 for Df = 2.6, a value remarkably similar to the 
Salpeter (1955) value, suggesting that the stellar initial mass function could be intimately 
related to the internal structure of molecular cloud complexes. 

In summary, we have derived the properties of fractal cloud complexes, no mattering 
the physical processes generating the fractal structure. It seems that Df ~ 2.6 is consistent 
with observations, but the relevance of this relatively high fractal dimension value has to be 
analyzed in future studies, mainly concerning the physical processes involved in the structure 
of the ISM. 
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